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Multifractal scaling analysis of nuclear giant resonance transition probability distri- 
butions is performed within the approximation which takes into account the one-particle- 
one-hole [Ip-lh) and 2p-2h states. A new measure to determine the fractal dimensions of 
the spectra is introduced. It is found that chaotic dynamics governing the decay leads to 
non-trivial multifractal scaling properties. Such a kind of scaling is absent in the case of 
regular dynamics. The degree of collectivity is another element which worsens the scaling. 
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1. Introduction 



There exist well established methods to identify chaotic behaviour in classical non- 
linear dynamical systems. It is, however, not fully clear what are the signatures of chaotic- 
ity at the quantum level. In fact, several indicators have been suggested. For systems 
with complex energy spectra the statistical analysis of their fluctuations has been intro- 
duced by Wigner in the context of compound nuclei. The Wigner form of the nearest 
neighbor spacing (NNS) distribution turns out to be a good indication that an underlying 
classical system is chaotic. This distribution can be simulated within the random matrix 
theory by the Gaussian orthogonal ensemble (GOE), while the regular systems' spectra 
have Poissonian distribution [1, 2]. The NNS test can be successfully applied to stationary 
systems and the information about the energy spectrum is sufficient. However, one would 
be interested to take into account more information included in the wave function [3] . The 
aim of this research is to find further signatures of chaos coming from the structure of 
the wave functions. In particular, chaotic dynamics on the classical level is known to be 
associated with various fractal structures obeying the scaling relations [4]. It would be 
extremely interesting to identify a trace of such structures also on the quantum level. 

To this end we shall consider the nuclear giant resonances as physically interesting 
excitations and we shall investigate their energy spectra and the transition probabilities re- 
sulting from the coupling to more complex configurations. These quantities are in principle 
measurable experimentally thus our predictions are verifiable. In particular, our calcula- 
tions have been done for the '^^Ca nucleus and its = 2+ excitations to have the well 
defined angular momentum and parity quantum numbers. The centroid energies of nuclear 
giant resonances are relatively well described within the mean field 1-particle-l-hole {Ip- 
Ih) approximation. To describe its decay one needs to include states of the np-nh type. 
Truncating on the 2p-2h level is, however, technically necessary (to have the hamiltonian 
matrix of a manageable size) and physically well justified [5]. As the nuclear forces are 
predominantly two-body in nature an initially excited giant resonance (a superposition of 
the Ip-lh states) can couple directly to the 2p-2h states only [6]. In particular, within 
this approximation one obtains a good agreement with experimental data. Furthermore, 
the corresponding spectral fiuctuations fulfill the GOE characteristics [7] which implies an 

2 



underlying chaotic dynamics and the fact that already the subspace of the 2p-2h states 
properly simulates relevant characteristics of the compound nuclei [8] . 



2. The physical model 



The general form of the Hamiltonian in our model reads: 



Vij,kl aja^j ai ak , 



(1) 



ij,kl 



where the first term (Hq) is the mean field part and the second term (V) is the residual 
interaction. For calculations we specify Hq in terms of a local Woods-Saxon potential 
including the Coulomb interaction. The interaction part, V, was taken as the zero range 
Landau-Migdal interaction with the empirical parameters taken from [9]. 

The Hilbert space in our model is spanned by the Ip-lh and 2p-2h vectors defined in 
terms of the creation and annihilation operators: 



|1) = alah\0) , |2) = al^al^tth^ahA^) , 



(2) 



such that they diagonalize the operator Hq in Ip-lh and 2p-2h sectors and (l|i7o|2) = 0. 
We prediagonalize the Hamiltonian H changing the basis: 



(3) 



and the Schrodinger equation then reads: 



Ei ^12 
^12 E^ 



= E 



(4) 



where Ei ^ are diagonal block-matrices, A12, Ai~ are the off-diagonal blocks and 



Ar, = J2cl{i\v\2)cl 



(5) 
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The A blocks are responsible for the coupling between these two sectors and 'leaking' of 
the probability from the Ip-lh sector to the 2p-2h one. 
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We consider two distinct and physically interesting cases: 

A) . No residual interaction in the subspace spanned by |2) (V = in that subspace). In 

this case the only non- vanishing matrix elements of the type {2\H\2') are the diagonal 
ones {{2\H\2') = {2\Ho\2') = €2522')- The fluctuation properties of 62 are then those 
characteristic for the regular systems [5]. 

B) . All the matrix elements are included which results in the GOE fluctuations of the 

energy spectra €2 typical for classically chaotic systems [5] . 
By comparing these two versions of our model one can then study the influence of chaotic 
dynamics on quantum decay. 

The transition probabilities that will be used in the following Section to deflne our 
measure (11) are taken as transitions between the ground state (|0)) and the excited states 

enumerated hy i = 1,2, . . . , N. In our notation a state \i) can be expressed as a 
superposition of states from both sectors: 

\i) = E«i Il) + E«2 |2) , (6) 
1 2 

and the quadruple transition operator reads: 

O = Q Y2{Q) , (7) 

where Q is the standard charge operator: Q = ^(l+ra), ts being the (isospin) Pauli matrix. 
The flrst term in Q deflnes the isoscalar and the second term the isovector component of 
the transition. 

The transition probabilities can be computed directly from the following formula: 

Pi = md\i)\'. (8) 

Since the transition operator deflned by eq. (7) is one-body in nature it picks-up the Ip-lh 
components of \i) only and thus: 

p, = I EXiCj(0|O|l)|'. (9) 
11 

In principle, any state \i) includes certain admixture of |1) and, therefore, an originally 
(with no coupling to |2)) localized transition strength becomes much more fragmented. 
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To have good quantum numbers we take the Ca nucleus with the angular momentum 
and parity: J"^ = 2+. Taking into account four mean field nuclear shells (two shells above 
and two below the Fermi surface) this implies 26 Ip-lh states and 3 014 2p-2h states. Here, 
the centroid energy is 31 and 25 MeV for isovector and isoscalar transitions respectively, 
with dispersion of the order of 5 MeV. Figs. 1 and 2 display the transition probabilities for 
different energy levels in the regular (A) and chaotic case (B) for isovector and isoscalar 
transitions, respectively. Even at this level one can observe a clear distinction between 
regular and chaotic cases. The later case even suggests a certain kind of self-similarity 
regarding the clustering and the relative size of the transitions. For this reason the following 
Section is an attempt to perform a more systematic analysis of these transition probabilities 
in the spirit of the Renyi exponents [10]. 



3. Scaling analysis 

The standard definition of a measure to compute the fractal (Renyi) dimensions (ex- 
ponents) for a set of n points in N boxes (intervals) of a size I is: 

71 ■ 

Pi{l) = - , (10) 
n 

where Ui = ni{l) denotes number of points in the i-th box (interval) and the measure is 
properly normalized due to: = n. With this measure one can determine the fractal 

dimensions dg (q being any real number) by taking the limit Z — > and extracting the 
exponent from the following formula [11]: 

J2pm ^ i^'^-'^'^ {1^0), (11) 

i 

where the log-log plot of the left hand side of eq. (11) vs. N is used to compute the 
standard scaling exponents dq. 

Applying these ideas to our energy levels only we get not a very interesting result as 
the energy spectrum of a typical quantum systems has dimension equal to one and the 
same holds true in both our cases. (There exist however systems whose energy spectra 
display a fractal character [12].) Therefore, consistently with our previous discussion we 
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are going to consider the sets of pairs of points, the i-th pair containing two numbers: the 
energy and the transition probabihty, {Ei,pi}. Taking into account the probabihties pi 
we have a more comphcated structure. Because to each energy corresponds exactly one 
probability the whole structure is expected to have dimension in the range [0, 1] if the 
idea of scaling applies. Of course, due to non-uniformity of the probability distributions 
the dimensions should depend on q [13]. 

The specific value of pi can be interpreted as a frequency with which the energy {Ei} 
is 'visited' and this modifies an effective number of energy levels in a given box. In this way 
each energy point gets a different weight. To have a measure with the proper normalization 
we define, instead of (10), the following new measure Pi{l): 



m) 



-\ -1 



J2 

L all 



E (12) 



X 

EiE i— th box 



where the summation in numerator goes over the probabilities whose energies are included 
in the i-th box. Here again the measure Pi{l) is properly normalized: X^j-Pi(0 = 1. The 
scaling exponent ("fractal dimension") with the measure (12) we will denote by Dg, while 
dq we reserve for the standard fractal dimension. From (12) it is clear that in the limit 
g — > (capacity dimension) the different nonzero probabilities pi give the same contribution 
to the left hand side of (11), because the measure Pi{l) is in the power of 0: 

limP.^fn = limp^f/) = I 1 z-th box is not empty, 

q^Q * ^ ' g^o * I if i-th box is empty. ^ ' 

Hence, we have: Dq = do, the last being the standard capacity dimension. In general, both 

measures and the corresponding dimensions differ for g > according to: < Dq < dq. 

The differences reflect a degree of non-uniformity in the probability distribution. To 

determine the scaling exponents we use the formula (11) with the new measure (12): 

Xq{l) ^ - (^^0)' (14) 

i 

and the log-log plot of the quantity XgiO of ^'i- (14) vs. the number of boxes N is used 
to extract the scaling exponents Dg in Figs. 3 and 4. 

The input data (of Figs. 1 and 2) consist of the order of ~ 2^^ data points, the 
number sufficient to display an exponential scaling but one should have in mind that some 
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statistical errors will be present, as the fractal dimension formula involves the I — > (or, 
equivalently, the N — > oo) limit. In fact, for the chaotic case we have got a fairly good 
scaling in the range of about 8 points in the log-log plot as can be seen from Figs. 3(B). 
What is important, this scaling persists up to the high q values, i.e. to the region where the 
structure of the probabilities is mostly probed. This kind of scaling has been considerably 
worsened in the regular case (Fig. 3(A) ). In the special case of the capacity dimension 
we get Do — 1, as in this case the scaling exponents are determined solely by the energy 
distribution (see eq. (2) ). This is not a very interesting limit and has not been plotted 
in Fig. 3(B). The regular case is also plotted for comparison, even though the scaling 
soon degrades with increasing q and choice of the scaling exponents is difficult in this case. 
Hence, for g > 2 the linear fits have not been plotted in Fig. 3(A). Also, as has been 
expected, we get larger differences between Dq and dg = 1 for greater g's. 

Analogous plots for the isoscalar transitions are displayed in Fig. 4. One can observe 
that in this case the scaling exponents (fractal dimensions) are slightly (about 5%) lower 
which is consistent with the less uniform distribution of the corresponding transition prob- 
abilities shown in Fig. 2. Also, the scaling is less evident (its range is shorter, see Fig. 4(B) ) 
even though the fluctuations in the subspace 2 remain the same. This may reflect the fact 
that the isoscalar giant resonance is more collective^ than its isovector counterpart and as 
such it is more resistive against decay [14]. In other words the collectivity is expected to 
regularize the dynamics and the result seems to confirm such a conjecture. 



4. Summary and conclusions 

In this paper we have investigated the scaling properties of nuclear giant quadrupole 
resonance transition probabilities in "^^Ca nucleus. These quantities have been computed 
within the approximation which includes the Ip-lh and 2p-2h states. The results still 
preserve the Wigner form of the NNS distribution of energy levels and they are in good 

^ The notion of the collectivity used here means localization of the transition strength 
in energy and should not be confused for the same notion used in the context of self- 
organization. 
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agreement with the experimental data (see Figs. 1, 2 and [5, 7]). 

To estimate the scahng properties we have introduced a measure defined by (12) which 
combines the energy spectra and the transition probabihties. This definition allows to make 
use of the concept of the generalized Renyi exponents [10, 11]. Based on this concept it has 
been shown that for both, the isovector and the isoscalar transitions, one can speak about 
the scaling exponents Dq of the multifractal type. This observation applies, however, only 
to the case when the physics of fragmentation is governed by the chaotic dynamics. These 
results can be treated as an interesting indication of what are the further signatures of 
classical chaos on the quantum level. 

For the isovector transitions the scaling is somewhat better and the exponents are 
about 5% higher. This effect may have to do with the fact that the isoscalar resonance 
is more collective than its isovector counterpart. The collectivity is a natural element 
regularizing the dynamics. 

It has been suggested [15] that dynamical systems of the type discussed in this paper 
can be simulated by the binary, self-similar and conservative random fragmentation process 
which yields universal behaviour independent of the precise fragmentation mechanism [15, 
16]. This gives another justification for neglecting the higher order excitations. 

We thank Marek Ploszajczak for very useful discussions. This research was supported 
by KBN Grant 2 P03B 140 10. 
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Figure captions. 



Fig. 1 Isovector quadrupole strength distribution in °Ca for the regular (A) and the chaotic 

(B) case. Note different energy scales. 
Fig. 2 The same as Fig. 1 for the isoscalar strength distribution. 

Fig. 3 The log-log plot of Xg(0 of ^Q- (14) vs. the number of boxes N (equivalent to the 
inverse of the box size I). Xgi^) is determined by the isovector quadruple transition 
probabilities in the regular (A) and chaotic (B) case for q = 2,4, 6. 

Fig. 4 The same as Fig. 3 for the isoscalar transition probabilities. 
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